home *** CD-ROM | disk | FTP | other *** search
/ The X-Philes (2nd Revision) / The X-Philes Number 1 (1995).iso / xphiles / hp48hor2 / tyko.src < prev    next >
Text File  |  1992-01-11  |  17KB  |  919 lines
  1. %%HP: T(3)A(D)F(.);
  2. @ TYKO by K. M. Sinenmaa
  3. DIR
  4.  MAIN
  5. \<< RCLF 'FL' STO DEG
  6. -42 SF MISC
  7.   IF 1 ==
  8.   THEN
  9. [ .606434 .374897 .259091 .827362 .779072 .993126 .700695 .485541
  10. .566441 .505498 .140023 .292498 .987353 .053856 .849694 .089608
  11. .056531 .814794 .133295 .882987 .821218 .870169 .400589 .664614
  12. .846912 .725368 .480856 .663854 .04102 .357355 ]
  13. [ .03660110129 .03629164709 .0367481952 .03386319198 .00273790931
  14. .0027377785 .011367714 .01136759566 .01136762384 .00445046867
  15. .00445036173 .00445040017 .00145575328 .00145561327 .00145569465
  16. .00023080893 .00023080893 .00023080893 .00009294371 .00009294371
  17. .00009294371 .00003269438 .00003269438 .00003265562 .00001672092
  18. .00001672092 .00001663715 .00001115482 .00001104864 .00001104864 ]
  19. C * + OBJ\-> DROP {
  20. F9 G9 L9 F8 G8 L8
  21. F7 G7 L7 F6 G6 L6
  22. F5 G5 L5 F4 G4 L4
  23. F2 G2 L2 F1 G1 L1
  24. GS LS D FM GM LM }
  25. 'BZ' STO 1 30
  26.     FOR j FP 360 *
  27. BZ j GET STO
  28.     NEXT CALS OBJ\->
  29. DROP EVAL 'r' STO
  30. EVAL 3600 / LS +
  31. F36 '\GlS' STO 1400
  32. .12 BEEP 1000 .2
  33. BEEP { { "Merc"
  34.     \<< CAL1 OBJ\->
  35. DROP EVAL 'R' STO
  36. EVAL '\Gb' STO EVAL
  37. 3600 / L1 + '\Gl' STO
  38. 2439 -.4 87 G1 * 1
  39. CONT
  40.     \>> } { "Venu"
  41.     \<< CAL2 OBJ\->
  42. DROP EVAL 'R' STO
  43. EVAL '\Gb' STO EVAL
  44. 3600 / L2 + '\Gl' STO
  45. 6050 -4 224 G2 * 2
  46. CONT
  47.     \>> } { "Mars"
  48.     \<< CAL4 OBJ\->
  49. DROP EVAL 'R' STO
  50. EVAL '\Gb' STO EVAL
  51. 3600 / L4 + '\Gl' STO
  52. 3397 -1.2 686 G4 *
  53. 3 CONT
  54.     \>> } { "Jupit"
  55.     \<< CAL5 OBJ\->
  56. DROP EVAL 'R' STO
  57. EVAL '\Gb' STO EVAL
  58. 3600 / L5 + '\Gl' STO
  59. 71398 -9.4 4360 G5
  60. * 4 CONT
  61.     \>> } { "Satur"
  62.     \<< CAL6 OBJ\->
  63. DROP EVAL 'R' STO
  64. EVAL '\Gb' STO EVAL
  65. 3600 / L6 + '\Gl' STO
  66. 60000 -9.4 10740 G6
  67. * 5 CONT
  68.     \>> } { "Uran"
  69.     \<< CAL7 OBJ\->
  70. DROP EVAL 'R' STO
  71. EVAL '\Gb' STO EVAL
  72. 3600 / L7 + '\Gl' STO
  73. 26320 -7.15 30690
  74. G7 * 6 CONT
  75.     \>> } { "Neptu"
  76.     \<< CAL8 OBJ\->
  77. DROP EVAL 'R' STO
  78. EVAL '\Gb' STO EVAL
  79. 3600 / L8 + '\Gl' STO
  80. 24300 -6.83 60181
  81. G8 * 7 CONT
  82.     \>> } { "Pluto"
  83.     \<< CAL9 OBJ\->
  84. DROP EVAL 'R' STO
  85. EVAL '\Gb' STO EVAL
  86. 3600 / L9 + '\Gl' STO
  87. 1500 -1 90480 G9 *
  88. 8 CONT
  89.     \>> } { "Moon"
  90.     \<< CALM OBJ\->
  91. DROP EVAL SIN INV
  92. DUP 'R' STO '\GD' STO
  93. EVAL 3600 / DUP '\Gb'
  94. STO 'b' STO EVAL
  95. 3600 / LM + F36 DUP
  96. '\Gl' STO 'l' STO
  97. 1738 27 GM * 9 CONT
  98.     \>> } { "Sun"
  99.     \<< r DUP 'R' STO
  100. '\GD' STO \GlS DUP '\Gl'
  101. STO 'l' STO 0 '\Gb'
  102. STO 0 'b' STO
  103. 696000 10 CONT
  104.     \>> } } TMENU
  105. "   Select  Object"
  106. PROMPT 'n' STO
  107. CLLCD
  108. "WAIT ..again" 1
  109. DISP
  110.     IF n 10 <
  111.     THEN LIST OBJ\->
  112. 7 DROPN SWAP 10000
  113. * + 1000000 / +
  114. SWAP 360 / NEG
  115. DATE+ '\Gt' STO
  116.     END
  117.     IF n 9 <
  118.     THEN \Gb 3600 /
  119. '\Gb' STO \Gl F36 '\Gl'
  120. STO GEOC 5 R \GD * 1
  121. l \Gl - COS + .5 * \v/
  122. / LOG * + 'm' STO
  123.     END
  124.     IF n 9 ==
  125.     THEN 2 * \GD 6378
  126. * / R\->D \->HMS
  127.     ELSE 2 * \GD
  128. 149597870 * / R\->D
  129. \->HMS
  130.     END 10000 *
  131.     IF n 8 >
  132.     THEN .6 *
  133.     END 'Q' STO
  134.   ELSE 1 'n' STO {
  135. { "STARS"
  136.     \<< NEWS
  137.     \>> } { "" } {
  138. "ESC"
  139.     \<< CONT
  140.     \>> } } TMENU
  141. "Browsing ..."
  142. PROMPT 1.396971 t *
  143. .0003086 t 2 ^ * +
  144. .013056 t *
  145. .0000092 t 2 ^ * -
  146. 5.12362 .241614 t *
  147. .0001122 t 2 ^ * +
  148. + DATA n GET OBJ\->
  149. DROP 'm' STO 'R'
  150. STO OVER 4 PICK +
  151. SIN 5 PICK * OVER +
  152. DUP '\Gb' STO 'b' STO
  153. TAN OVER 4 ROLL +
  154. COS * ROT * NEG ROT
  155. + + DUP '\Gl' STO 'l'
  156. STO n 10 + 'n' STO
  157.   END CONV AAUT { Y
  158. M LT \GlM \O/ TZ C T t
  159. \Ge LM GM FM D \GWM LS
  160. GS L1 G1 F1 L2 G2 F2
  161. L4 G4 F4 L5 G5 F5
  162. L6 G6 F6 L7 G7 F7
  163. L8 G8 F8 L9 G9 F9 T
  164. BZ \GlS h l b } PURGE
  165. ERASE NM n GET 2
  166. \->GROB PICT STO PICT
  167. RCL { # 0d # 58d }
  168. LIST K ZZ 6 FIX {
  169. # 60d # 8d }
  170. "GST 0UT= " UT \->HMS
  171. ZZ { # 60d # 15d }
  172. "LST= " ST \->HMS ZZ
  173. { # 0d # 10d }
  174. "\Gl= " \Gl \->HMS ZZ {
  175. # 0d # 16d } "\Gb= "
  176. \Gb \->HMS ZZ { # 0d
  177. # 24d } "r= " R ZZ
  178. { # 0d # 30d }
  179. "\GD= " \GD ZZ { # 0d
  180. # 38d } "\Gd= " \Gd
  181. \->HMS ZZ { # 0d
  182. # 44d } "\Ga= " \Ga 15
  183. / \->HMS ZZ { # 60d
  184. # 23d } "AZ= " A
  185. \->HMS ZZ { # 60d
  186. # 29d } "a = " a
  187. \->HMS ZZ { # 0d
  188. # 51d } "\Gt= " \Gt ZZ
  189. { # 60d # 44d }
  190. "LT\|^= " RI ZZ {
  191. # 60d # 51d }
  192. "LT\|v= " SE ZZ 2 FIX
  193. { # 57d # 37d }
  194. "D= " Q ZZ { # 98d
  195. # 37d } "m= " m ZZ
  196. { # 68d # 1d }
  197. "JDAY= " JD ZZ PICT
  198. STO 400 .2 BEEP 800
  199. .2 BEEP 1600 .3
  200. BEEP PICT RCL \->LCD
  201. { } PVIEW FL STOF {
  202. n \GD \Gd \Ga \Gl \Gb r R FL
  203. A a ST \Gt m Q UT JD
  204. RI SE } PURGE 2
  205. MENU
  206. \>>
  207. MISC
  208. \<< POP STD { {
  209. "DATE"
  210.   \<<
  211. "DATE-format DD.MMYYYY"
  212. ":DATE:" INPUT OBJ\->
  213. DUP IP SWAP FP 100
  214. * DUP IP SWAP FP
  215. 10000 * LIST { 1 3
  216. } ROT PUT { 1 2 }
  217. ROT PUT { 1 1 } ROT
  218. PUT 'LIST' STO NEW
  219.   \>> } { "TIME"
  220.   \<<
  221. "TIME-format HH.MMSS"
  222. ":TIME:"
  223. INPUT OBJ\-> LIST { 3
  224. 1 } ROT PUT 'LIST'
  225. STO NEW
  226.   \>> } { "LONG"
  227.   \<<
  228. "Longitude in decimals"
  229. "
  230. west=neg
  231. " +
  232. ":\Gl:" INPUT OBJ\->
  233. LIST { 2 1 } ROT
  234. PUT 'LIST' STO NEW
  235.   \>> } { "LATIT"
  236.   \<<
  237. "Latitude in decimals"
  238. ":\O/:" INPUT OBJ\->
  239. LIST { 2 2 } ROT
  240. PUT 'LIST' STO NEW
  241.   \>> } { "ZONE"
  242.   \<<
  243. "Integer -12 to 12"
  244. ":Time Zone:" INPUT
  245. OBJ\-> LIST { 2 3 }
  246. ROT PUT 'LIST' STO
  247. NEW
  248.   \>> } { "ESC"
  249.   \<< { { "SUN.S"
  250.     \<< 1 CONT
  251.     \>> } { "" } {
  252. "STARS"
  253.     \<< 2 CONT
  254.     \>> } } TMENU
  255.   \>> } } TMENU NEW
  256. HALT CLLCD
  257. "WAIT ..." 1 DISP
  258. LIST OBJ\-> 3 DROPN
  259. 'LT' STO 'TZ' STO
  260. '\O/' STO '\GlM' STO
  261. 'Y' STO 'M' STO 367
  262. Y * 7 Y M 9 + 12 /
  263. IP + * 4 / IP - 3 Y
  264. M 9 - 7 / IP + 100
  265. / IP 1 + * 4 / IP -
  266. 275 M * 9 / IP +
  267. 1721028.5 + + DUP
  268. LT TZ - F24 HMS\-> 24
  269. / + 'JD' STO
  270. 2451545 - 36525 /
  271. 8640184.81287 *
  272. 24110.54841 + JD
  273. 2451545 - DUP 'C'
  274. STO 36525 / 't' STO
  275. .093104 t 2 ^ * +
  276. .0000062 t 3 ^ * -
  277. 3600 / F24 DUP 'UT'
  278. STO 1.00273790935
  279. LT TZ - F24 HMS\-> *
  280. + \GlM 15 / + F24
  281. 'ST' STO t 1 + 'T'
  282. STO .347343 C
  283. .00014709391 * -
  284. '\GWM' STO 84428 47 T
  285. * - 9 \GWM COS * +
  286. 3600 / '\Ge' STO
  287. \>>
  288. CALM
  289. \<< POP { '22640*SIN(
  290. GM)-4586*SIN(GM-2*D
  291. )+2370*SIN(2*D)+769
  292. *SIN(2*GM)-668*SIN(
  293. GS)-412*SIN(2*FM)-
  294. 212*SIN(2*GM-2*D)-
  295. 206*SIN(GM-2*D+GS)+
  296. 192*SIN(GM+2*D)+165
  297. *SIN(2*D-GS)+148*
  298. SIN(GM-GS)-125*SIN(
  299. D)-110*SIN(G2+GS)-
  300. 55*SIN(2*FM-2*D)-45
  301. *SIN(GM+2*FM)+40*
  302. SIN(GM-2*FM)-38*SIN
  303. (GM-4*D)+36*SIN(3*
  304. GM)-31*SIN(2*GM-4*D
  305. )+28*SIN(GM-2*D-GS)
  306. -24*SIN(2*D+GS)+19*
  307. SIN(GM-D)+18*SIN(D+
  308. GS)+15*SIN(GM+2*D-
  309. GS)+14*SIN(2*GM+2*D
  310. )+14*SIN(4*D)-13*
  311. SIN(3*GM-2*D)-11*
  312. SIN(GM+16*LS-18*L2)
  313. +10*SIN(2*GM-GS)+9*
  314. SIN(GM-2*FM-2*D)+9*
  315. COS(GM+16*LS-18*L2)
  316. -9*SIN(2*GM-2*D+GS)
  317. -8*SIN(GM+D)+8*SIN(
  318. 2*D-2*GS)-8*SIN(2*
  319. GM+GS)-7*SIN(2*GS)-
  320. 7*SIN(GM-2*D+2*GS)+
  321. 7*SIN(\GWM)-6*SIN(GM-
  322. 2*FM+2*D)-6*SIN(2*
  323. FM+2*D)-4*SIN(GM-4*
  324. D+GS)+4*T*COS(GM+16
  325. *LS-18*L2)-4*SIN(2*
  326. GM+2*FM)+4*T*SIN(GM
  327. +16*LS-18*L2)+3*SIN
  328. (GM-3*D)-3*SIN(GM+2
  329. *D+GS)-3*SIN(2*GM-4
  330. *D+GS)+3*SIN(GM-2*
  331. GS)+3*SIN(GM-2*D-2*
  332. GS)-2*SIN(2*GM-2*D-
  333. GS)-2*SIN(2*FM-2*D+
  334. GS)+2*SIN(GM+4*D)+2
  335. *SIN(4*GM)+2*SIN(4*
  336. D-GS)+2*SIN(2*GM-D)
  337. ' '18461*SIN(FM)+
  338. 1010*SIN(GM+FM)+
  339. 1000*SIN(GM-FM)-624
  340. *SIN(FM-2*D)-199*
  341. SIN(GM-FM-2*D)-167*
  342. SIN(GM+FM-2*D)+117*
  343. SIN(FM+2*D)+62*SIN(
  344. 2*GM+FM)+33*SIN(GM-
  345. FM+2*D)+32*SIN(2*GM
  346. -FM)-30*SIN(FM-2*D+
  347. GS)-16*SIN(2*GM+FM-
  348. 2*D)+15*SIN(GM+FM+2
  349. *D)+12*SIN(FM-2*D-
  350. GS)-9*SIN(GM-FM-2*D
  351. +GS)-8*SIN(FM+\GWM)+8
  352. *SIN(FM+2*D-GS)-7*
  353. SIN(GM+FM-2*D+GS)+7
  354. *SIN(GM+FM-GS)-7*
  355. SIN(GM+FM-4*D)-6*
  356. SIN(FM+GS)-6*SIN(3*
  357. FM)+6*SIN(GM-FM-GS)
  358. -5*SIN(FM+D)-5*SIN(
  359. GM+FM+GS)-5*SIN(GM-
  360. FM+GS)+5*SIN(FM-GS)
  361. +5*SIN(FM-D)+4*SIN(
  362. 3*GM+FM)-4*SIN(FM-4
  363. *D)-3*SIN(GM-FM-4*D
  364. )+3*SIN(GM-3*FM)-2*
  365. SIN(2*GM-FM-4*D)-2*
  366. SIN(3*FM-2*D)+2*SIN
  367. (2*GM-FM+2*D)+2*SIN
  368. (GM-FM+2*D-GS)+2*
  369. SIN(2*GM-FM-2*D)+2*
  370. SIN(3*GM-FM)' '
  371. .9508+.0518*COS(
  372. 134.9+447198.85*T)+
  373. .0095*COS(259.2-
  374. 413335.38*T)+.0078*
  375. COS(235.7+890534.23
  376. *T)+.0028*COS(269.9
  377. +954397.7*T)' }
  378. \>>
  379. CAL1
  380. \<< POP { '84378*SIN(
  381. G1)+10733*SIN(2*G1)
  382. +1892*SIN(3*G1)-646
  383. *SIN(2*F1)+381*SIN(
  384. 4*G1)-306*SIN(G1-2*
  385. F1)-274*SIN(G1+2*F1
  386. )-92*SIN(2*G1+2*F1)
  387. +83*SIN(5*G1)-28*
  388. SIN(3*G1+2*F1)+25*
  389. SIN(2*G1-2*F1)+19*
  390. SIN(6*G1)-9*SIN(4*
  391. G1+2*F1)+8*T*SIN(G1
  392. )+7*COS(2*G1-5*G2)' 
  393. '24134*SIN(F1)+5180
  394. *SIN(G1-F1)+4910*
  395. SIN(G1+F1)+1124*SIN
  396. (2*G1+F1)+271*SIN(3
  397. *G1+F1)+132*SIN(2*
  398. G1-F1)+67*SIN(4*G1+
  399. F1)+18*SIN(3*G1-F1)
  400. +17*SIN(5*G1+F1)-10
  401. *SIN(3*F1)-9*SIN(G1
  402. -3*F1)' '.39528-
  403. .07834*COS(G1)-
  404. .00795*COS(2*G1)-
  405. .00121*COS(3*G1)-
  406. .00022*COS(4*G1)' }
  407. \>>
  408. CAL2
  409. \<< POP { '2814*SIN(
  410. G2)-181*SIN(2*F2)-
  411. 20*T*SIN(G2)+12*SIN
  412. (2*G2)-10*COS(2*GS-
  413. 2*G2)+7*COS(3*GS-3*
  414. G2)' '12215*SIN(F2)
  415. +83*SIN(G2+F2)+83*
  416. SIN(G2-F2)' '.72335
  417. -.00493*COS(G2)' }
  418. \>>
  419. CAL4
  420. \<< POP { '38451*SIN(
  421. G4)+2238*SIN(2*G4)+
  422. 181*SIN(3*G4)-52*
  423. SIN(2*F4)+37*T*SIN(
  424. G4)-22*COS(G4-2*G5)
  425. -19*SIN(G4-G5)+17*
  426. COS(G4-G5)+17*SIN(4
  427. *G4)-16*COS(2*G4-2*
  428. G5)+13*COS(GS-2*G4)
  429. -10*SIN(G4-2*F4)-10
  430. *SIN(G4+2*F4)+7*COS
  431. (GS-G4)-7*COS(2*GS-
  432. 3*G4)-5*SIN(G2-3*G4
  433. )-5*SIN(GS-G4)-5*
  434. SIN(GS-2*G4)-4*COS(
  435. 2*GS-4*G4)+4*T*SIN(
  436. 2*G4)+4*COS(G5)+3*
  437. COS(G2-3*G4)+3*SIN(
  438. 2*G4-2*G5)' '6603*
  439. SIN(F4)+622*SIN(G4-
  440. F4)+615*SIN(G4+F4)+
  441. 64*SIN(2*G4+F4)' '
  442. 1.53031-.1417*COS(
  443. G4)-.0066*COS(2*G4)
  444. -.00047*COS(3*G4)'
  445. }
  446. \>>
  447. CAL5
  448. \<< POP { '19934*SIN(
  449. G5)+5023*T+2511+
  450. 1093*COS(2*G5-5*G6)
  451. +601*SIN(2*G5)-479*
  452. SIN(2*G5-5*G6)-185*
  453. SIN(2*G5-2*G6)+137*
  454. SIN(3*G5-5*G6)-131*
  455. SIN(G5-2*G6)+79*COS
  456. (G5-G6)-76*COS(2*G5
  457. -2*G6)-74*T*COS(G5)
  458. +68*T*SIN(G5)+66*
  459. COS(2*G5-3*G6)+63*
  460. COS(3*G5-5*G6)+53*
  461. COS(G5-5*G6)+49*SIN
  462. (2*G5-3*G6)-43*T*
  463. SIN(2*G5-5*G6)-37*
  464. COS(G5)+25*SIN(2*L5
  465. )+25*SIN(3*G5)-23*
  466. SIN(G5-5*G6)-19*T*
  467. COS(2*G5-5*G6)+17*
  468. COS(2*G5-4*G6)+17*
  469. COS(3*G5-3*G6)-14*
  470. SIN(G5-G6)-13*SIN(3
  471. *G5-4*G6)-9*COS(2*
  472. L5)+9*COS(G6)-9*SIN
  473. (G6)-9*SIN(3*G5-2*
  474. G6)+9*SIN(4*G5-5*G6
  475. )+9*SIN(2*G5-6*G6+3
  476. *G7)-8*COS(4*G5-10*
  477. G6)+7*COS(3*G5-4*G6
  478. )-7*COS(G5-3*G6)-7*
  479. SIN(4*G5-10*G6)-7*
  480. SIN(G5-3*G6)+6*COS(
  481. 4*G5-5*G6)-6*SIN(3*
  482. G5-3*G6)+5*COS(2*G6
  483. )-4*SIN(4*G5-4*G6)-
  484. 4*COS(3*G6)+4*COS(2
  485. *G5-G6)-4*COS(3*G5-
  486. 2*G6)-4*T*COS(2*G5)
  487. +3*T*SIN(2*G5)+3*
  488. COS(5*G6)+3*COS(5*
  489. G5-10*G6)+3*SIN(2*
  490. G6)-2*SIN(2*L5-G5)+
  491. 2*SIN(2*L5+G5)-2*T*
  492. SIN(3*G5-5*G6)-2*T*
  493. SIN(G5-5*G6)' '
  494. -4692*COS(G5)+259*
  495. SIN(G5)+227-227*COS
  496. (2*G5)+30*T*SIN(G5)
  497. +21*T*COS(G5)+16*
  498. SIN(3*G5-5*G6)-13*
  499. SIN(G5-5*G6)-12*COS
  500. (3*G5)+12*SIN(2*G5)
  501. +7*COS(3*G5-5*G6)-5
  502. *COS(G5-5*G6)' '
  503. 5.20883-.25122*COS(
  504. G5)-.00604*COS(2*G5
  505. )+.0026*COS(2*G5-2*
  506. G6)-.0017*COS(3*G5-
  507. 5*G6)-.00106*SIN(2*
  508. G5-2*G6)-.00091*T*
  509. SIN(G5)-.00084*T*
  510. COS(G5)+.00069*SIN(
  511. 2*G5-3*G6)-.00067*
  512. SIN(G5-5*G6)+.00066
  513. *SIN(3*G5-5*G6)+
  514. .00063*SIN(G5-G6)-
  515. .00051*COS(2*G5-3*
  516. G6)-.00046*SIN(G5)-
  517. .00029*COS(G5-5*G6)
  518. +.00027*COS(G5-2*G6
  519. )-.00022*COS(3*G5)-
  520. .00021*SIN(2*G5-5*
  521. G6)' }
  522. \>>
  523. CAL6
  524. \<< POP { '23045*SIN(
  525. G6)+5014*T-2689*COS
  526. (2*G5-5*G6)+2507+
  527. 1177*SIN(2*G5-5*G6)
  528. -826*COS(2*G5-4*G6)
  529. +802*SIN(2*G6)+425*
  530. SIN(G5-2*G6)-229*T*
  531. COS(G6)-153*COS(2*
  532. G5-6*G6)-142*T*SIN(
  533. G6)-114*COS(G6)+101
  534. *T*SIN(2*G5-5*G6)-
  535. 70*COS(2*L6)+67*SIN
  536. (2*L6)+66*SIN(2*G5-
  537. 6*G6)+60*T*COS(2*G5
  538. -5*G6)+41*SIN(G5-3*
  539. G6)+39*SIN(3*G6)+31
  540. *SIN(G5-G6)+31*SIN(
  541. 2*G5-2*G6)-29*COS(2
  542. *G5-3*G6)-28*SIN(2*
  543. G5-6*G6+3*G7)+28*
  544. COS(G5-3*G6)+22*T*
  545. SIN(2*G5-4*G6)-22*
  546. SIN(G6-3*G7)+20*SIN
  547. (2*G5-3*G6)+20*COS(
  548. 4*G5-10*G6)+19*COS(
  549. 2*G6-3*G7)+19*SIN(4
  550. *G5-10*G6)-17*T*COS
  551. (2*G6)-16*COS(G6-3*
  552. G7)-12*SIN(2*G5-4*
  553. G6)+12*COS(G5)-12*
  554. SIN(2*G6-2*G7)-11*T
  555. *SIN(2*G6)-11*COS(2
  556. *G5-7*G6)+10*SIN(2*
  557. G6-3*G7)+10*COS(2*
  558. G5-2*G6)+9*SIN(4*G5
  559. -9*G6)-8*SIN(G6-2*
  560. G7)-8*COS(2*L6+G6)+
  561. 8*COS(2*L6-G6)+8*
  562. COS(G6-G7)-8*SIN(2*
  563. L6-G6)+7*SIN(2*L6+
  564. G6)-7*COS(G5-2*G6)-
  565. 7*COS(2*G6)-6*T*SIN
  566. (4*G5-10*G6)+6*T*
  567. COS(4*G5-10*G6)+6*T
  568. *SIN(2*G5-6*G7)-5*
  569. SIN(3*G5-7*G6)-5*
  570. COS(3*G5-3*G6)-5*
  571. COS(2*G6-2*G7)+5*
  572. SIN(3*G5-4*G6)+5*
  573. SIN(2*G5-7*G6)+4*
  574. SIN(3*G5-3*G6)+4*
  575. SIN(3*G5-5*G6)+4*T*
  576. COS(G5-2*G6)+3*T*
  577. COS(2*G5-4*G6)+3*
  578. COS(2*G5-6*G6+3*G7)
  579. -3*T*SIN(2*L6)+3*T*
  580. COS(2*G5-6*G6)-3*T*
  581. COS(2*L6)+3*COS(3*
  582. G5-7*G6)+3*COS(4*G5
  583. -9*G6)+3*SIN(3*G5-6
  584. *G6)+3*SIN(2*G5-G6)
  585. +3*SIN(G5-4*G6)+2*
  586. COS(3*G6-3*G7)+2*T*
  587. SIN(G5-2*G6)+2*SIN(
  588. 4*G6)-2*COS(3*G5-4*
  589. G6)-2*COS(2*G5-G6)-
  590. 2*SIN(2*G5-7*G6)+2*
  591. COS(G5-4*G6)+2*COS(
  592. 4*G5-11*G6)-2*SIN(
  593. G6-G7)' '8297*SIN(
  594. G6)-3346*COS(G6)+
  595. 462*SIN(2*G6)-189*
  596. COS(2*G6)+185+79*T*
  597. COS(G6)-71*COS(2*G5
  598. -4*G6)+46*SIN(2*G5-
  599. 6*G6)-45*COS(2*G5-6
  600. *G6)+29*SIN(3*G6)-
  601. 20*COS(2*G5-3*G6)+
  602. 18*T*SIN(G6)-14*COS
  603. (2*G5-5*G6)-11*COS(
  604. 3*G6)-10*T+9*SIN(G5
  605. -3*G6)+8*SIN(G5-G6)
  606. -6*SIN(2*G5-3*G6)+5
  607. *SIN(2*G5-7*G6)-5*
  608. COS(2*G5-7*G6)+4*
  609. SIN(2*G5-5*G6)-4*T*
  610. SIN(2*G6)-3*COS(G5-
  611. G6)+3*COS(G5-3*G6)+
  612. 3*T*SIN(2*G5-4*G6)+
  613. 3*SIN(G5-2*G6)+2*
  614. SIN(4*G6)-2*COS(2*
  615. G5-2*G6)' '9.55774-
  616. .53252*COS(G6)-
  617. .01878*SIN(2*G5-4*
  618. G6)-.01482*COS(2*G6
  619. )+.00817*SIN(G5-G6)
  620. -.00539*COS(G5-2*G6
  621. )-.00524*T*SIN(G6)+
  622. .00349*SIN(2*G5-5*
  623. G6)+.00347*SIN(2*G5
  624. -6*G6)+.00328*T*COS
  625. (G6)-.00225*SIN(G6)
  626. +.00149*COS(2*G5-6*
  627. G6)-.00126*COS(2*G5
  628. -2*G6)+.00104*COS(
  629. G5-G6)+.00101*COS(2
  630. *G5-5*G6)+.00098*
  631. COS(G5-3*G6)-.00073
  632. *COS(2*G5-3*G6)-
  633. .00062*COS(3*G6)+
  634. .00042*SIN(2*G6-3*
  635. G7)+.00041*SIN(2*G5
  636. -2*G6)-.0004*SIN(G5
  637. -3*G6)+.0004*COS(2*
  638. G5-4*G6)-.00028*T-
  639. .00023*SIN(G5)+
  640. .0002*SIN(2*G5-7*G6
  641. )' }
  642. \>>
  643. CAL7
  644. \<< POP { '19397*SIN(
  645. G7)+570*SIN(2*G7)-
  646. 536*T*COS(G7)+143*
  647. SIN(G6-2*G7)+110*T*
  648. SIN(G7)+102*SIN(G6-
  649. 3*G7)+76*COS(G6-3*
  650. G7)-49*SIN(G5-G7)+
  651. 64*T-30*T*COS(2*G7)
  652. +29*SIN(2*G5-6*G6+3
  653. *G7)+29*COS(2*G7-2*
  654. G8)-28*COS(G7-G8)+
  655. 23*SIN(3*G7)-21*COS
  656. (G5-G7)+20*SIN(G7-
  657. G8)+20*COS(G6-2*G7)
  658. -19*COS(G6-G7)+17*
  659. SIN(2*G7-3*G8)+14*
  660. SIN(3*G7-3*G8)+13*
  661. SIN(G6-G7)-24*T*COS
  662. (G7)-12*COS(G7)+10*
  663. SIN(2*G7-2*G8)-9*
  664. SIN(2*F7)-18*SIN(G7
  665. )+9*COS(2*G7-3*G8)+
  666. 8*T*COS(G6-2*G7)+7*
  667. T*COS(G6-3*G7)-7*T*
  668. SIN(G6-3*G7)+7*T*
  669. SIN(2*G7)+6*SIN(2*
  670. G5-6*G6+2*G7)+6*COS
  671. (2*G5-6*G6+2*G7)+5*
  672. SIN(G6-4*G7)-4*SIN(
  673. 3*G7-4*G8)+4*COS(3*
  674. G7-3*G8)-3*COS(G8)-
  675. 2*SIN(G8)' '2775*
  676. SIN(F7)+131*SIN(G7-
  677. F7)+130*SIN(G7+F7)' 
  678. '19.21216-.90154*
  679. COS(G7)-.02488*T*
  680. SIN(G7)-.02121*COS(
  681. 2*G7)-.00585*COS(G6
  682. -2*G7)-.00508*T*COS
  683. (G7)-.00451*COS(G5-
  684. G7)+.00336*SIN(G6-
  685. G7)+.00198*SIN(G5-
  686. G7)+.00118*COS(G6-3
  687. *G7)+.00107*SIN(G6-
  688. 2*G7)-.00103*T*SIN(
  689. 2*G7)-.00081*COS(3*
  690. G7-3*G8)' }
  691. \>>
  692. CAL8
  693. \<< POP { '3523*SIN(
  694. G8)-50*SIN(2*F8)-43
  695. *T*COS(G8)+29*SIN(
  696. G5-G8)+19*SIN(2*G8)
  697. -18*COS(G5-G8)+13*
  698. COS(G6-G8)+13*SIN(
  699. G6-G8)-9*SIN(2*G7-3
  700. *G8)+9*COS(2*G7-2*
  701. G8)-5*COS(2*G7-3*G8
  702. )-4*T*SIN(G8)+4*COS
  703. (G7-2*G8)+8*T*SIN(
  704. G8)' '6404*SIN(F8)+
  705. 55*SIN(G8+F8)+55*
  706. SIN(G8-F8)-33*T*SIN
  707. (F8)' '30.07175-
  708. .25701*COS(G8)-
  709. .00787*COS(2*L7-G7-
  710. 2*L8)+.00409*COS(G5
  711. -G8)-.00314*T*SIN(
  712. G8)+.0025*SIN(G5-G8
  713. )-.00194*SIN(G6-G8)
  714. +.00185*COS(G6-G8)'
  715. }
  716. \>>
  717. CAL9
  718. \<< POP { '101577*SIN
  719. (G9)+15517*SIN(2*G9
  720. )-3593*SIN(2*F9)+
  721. 3414*SIN(3*G9)-2201
  722. *SIN(G9-2*F9)-1871*
  723. SIN(G9+2*F9)+839*
  724. SIN(4*G9)-757*SIN(2
  725. *G9+2*F9)-285*SIN(3
  726. *G9+2*F9)+454*T*SIN
  727. (G9)+218*SIN(2*G9-2
  728. *F9)+200*T*SIN(G9)' 
  729. '57726*SIN(F9)+
  730. 15257*SIN(G9-F9)+
  731. 14102*SIN(G9+F9)+
  732. 3870*SIN(2*G9+F9)+
  733. 1138*SIN(3*G9+F9)+
  734. 472*SIN(2*G9-F9)+
  735. 353*SIN(4*G9+F9)-
  736. 144*SIN(G9-3*F9)-
  737. 119*SIN(3*F9)-111*
  738. SIN(G9+3*F9)' '
  739. 40.74638-9.58235*
  740. COS(G9)-1.16703*COS
  741. (2*G9)-.22649*COS(3
  742. *G9)-.04996*COS(4*
  743. G9)' }
  744. \>>
  745. CALS
  746. \<< POP { '6910*SIN(
  747. GS)+72*SIN(2*GS)-17
  748. *T*SIN(GS)-7*COS(GS
  749. -G5)+6*SIN(LM-LS)+5
  750. *SIN(4*GS-8*G4+3*G5
  751. )-5*COS(2*GS-2*G2)-
  752. 4*SIN(GS-G2)+4*COS(
  753. 4*GS-8*G4+3*G5)+3*
  754. SIN(2*GS-2*G2)-3*
  755. SIN(G5)-3*SIN(2*GS-
  756. 2*G5)' '1.00014-
  757. .01675*COS(GS)-
  758. .00014*COS(2*GS)' }
  759. \>>
  760. AAUT
  761. \<< POP -16 CF ST 15
  762. * \Ga - 'h' STO h COS
  763. \Gd COS \O/ COS * * \Gd
  764. SIN \O/ SIN * + ASIN
  765. 'a' STO h COS \Gd COS
  766. \O/ SIN * * \Gd SIN \O/
  767. COS * - h SIN \Gd COS
  768. * \->V2 -16 SF V\-> F36
  769. 'A' STO DROP \Gd TAN
  770. \O/ TAN * DUP
  771.   IF ABS 1 >
  772.   THEN DROP
  773.     IF a 0 >
  774.     THEN
  775. "Continuous"
  776.     ELSE "    No"
  777.     END
  778. "phenomenon"
  779.   ELSE NEG ACOS \Ga
  780. OVER - 15 / UT -
  781. SWAP \Ga + 15 / UT -
  782.     IF n 9 \=/
  783.     THEN \Gd ABS
  784.       IF \O/ 0 >
  785.       THEN \O/ + 9.4
  786. / SIN DUP ROT + ROT
  787. ROT -
  788.       ELSE \O/ NEG +
  789. 9.3 / SIN DUP ROT +
  790. ROT ROT -
  791.       END
  792.     ELSE \Gd \O/ ABS +
  793. SIN 2 / DUP ROT +
  794. ROT ROT +
  795.     END TZ 15 * \GlM
  796. - 15 / DUP ROT +
  797. F24 \->HMS SWAP ROT +
  798. F24 \->HMS
  799.   END 'SE' STO 'RI'
  800. STO
  801. \>>
  802. CONV
  803. \<< POP -16 CF l COS
  804. l SIN \Ge COS * b TAN
  805. \Ge SIN * - \->V2 -16
  806. SF V\-> F36 '\Ga' STO
  807. DROP b SIN \Ge COS *
  808. b COS \Ge SIN l SIN *
  809. * + ASIN '\Gd' STO
  810. \>>
  811. DATA
  812. { { 69.78901342
  813. -5.46743986 21 .85
  814. } { 301.75502956
  815. 29.30364784 5.1 .77
  816. } { 249.7616937
  817. -4.56997796 130 .96
  818. } { 204.2332694
  819. 30.7362084 11 -.04
  820. } { 88.838405
  821. -16.0275303 200 .01
  822. } { 81.85774154
  823. 22.86434254 14 .08
  824. } { 335.3285854
  825. 59.90627092 560
  826. 1.25 } {
  827. 115.78545193
  828. -16.019622324 3.5
  829. .38 } {
  830. 149.82870578
  831. .46469195 26 1.35 }
  832. { 76.82923859
  833. -31.12276919 280
  834. .12 } { 104.081575
  835. -39.6052516 2.7
  836. -1.46 } }
  837. GEOC
  838. \<< POP -16 CF R 2 ^
  839. r 2 ^ + 2 R r \GlS \Gl
  840. - 180 + COS * * * -
  841. \v/ '\GD' STO R \Gb SIN *
  842. \GD / ASIN 'b' STO \Gb
  843. COS \Gl COS R * * r
  844. \GlS COS * + R \Gb COS
  845. \Gl SIN * * r \GlS SIN
  846. * + \->V2 -16 SF V\->
  847. F36 'l' STO DROP
  848. \>>
  849. NEW
  850. \<<
  851. "      Changeable
  852. "
  853. "
  854. " + LIST \->STR +
  855. "
  856. " +
  857. "
  858.     ESC to continue
  859. "
  860. + CLLCD 1 DISP 3
  861. FREEZE
  862. \>>
  863. NEWS
  864. \<< " Number ?"
  865. ":No:" INPUT OBJ\->
  866. 'n' STO
  867.   IF n 11 >
  868.   THEN CLLCD
  869. "Not computationable"
  870. 1 DISP 3 FREEZE
  871.   ELSE CLLCD NM n
  872. 10 + GET 1 DISP 3
  873. FREEZE
  874.   END
  875. \>>
  876. NM
  877. { "Mercury" "Venus"
  878. "Mars" "Jupiter"
  879. "Saturn" "Uranus"
  880. "Neptune" "Pluto"
  881. "Moon" "Sun"
  882. "Aldebaran"
  883. "Altair" "Antares"
  884. "Arcturus"
  885. "Betelgeuse"
  886. "Capella" "Deneb"
  887. "Procyon" "Regulus"
  888. "Rigel" "Sirius" }
  889. LIST
  890. [[ 10 12 1991 ]
  891.  [ 24.954 60.162 2 ]
  892.  [ 2 0 0 ]]
  893. POP
  894. \<<
  895.   IF 'FL' VTYPE 0 <
  896.   THEN
  897. "Something lost ?"
  898. DOERR
  899.   END
  900. \>>
  901. ZZ
  902. \<< \->STR + 1 \->GROB
  903. GOR
  904. \>>
  905. F36
  906. \<< DUP 360 / IP
  907. 360 * -
  908.   IF DUP 0 <
  909.   THEN 360 +
  910.   END
  911. \>>
  912. F24
  913. \<< DUP 24 / IP 24 *
  914. -
  915.   IF DUP 0 <
  916.   THEN 24 +
  917.   END
  918. \>>
  919.